Optical control of entangled states in quantum wells 
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We present theory and calculations for coherent high-fidelity quantum control of many-particle 
states in semiconductor quantum wells. We show that coupling a two-electron double quantum dot 
to a terahertz optical source enables targeted excitations that are one to two orders of magnitude 
faster and significantly more accurate than those obtained with electric gates. The optical fields 
subject to physical constraints are obtained through quantum optimal control theory that we apply 
in conjunction with the numerically exact solution of the time-dependent Schrodinger equation. Our 
ability to coherently control arbitrary two-electron states, and to maximize the entanglement, opens 
up further perspectives in solid-state quantum information. 
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PACS numbers: 78.67.Hc, 73.21. La, 78.20.Bh, 03.67.Bg 

Quest for solid-state quantum computing^ in low- 
dimensional semiconductor materials has led to exper- 
imental breakthroughs in the initialization, processing, 
and readout of single- and two-particle states in cou- 
pled quantum dots 4^— The majority of proposed quan- 
tum computation schemes focus on spin control^ due to 
the long decoherence times up to microseconds. However, 
rapid development of light sources and ultrafast science in 
atomic physics 8 is likely to bring methods and techniques 
to coherent control of electric charge^ in low-dimensional 
systems. 

Electron dynamics in double quantum dots (DQDs) 
has been subject to various theoretical approaches. Be- 
yond few-level models, realistic DQDs have been con- 
trolled with gate voltage o 10 ' 11 ' 13 and optimized laser 
pulses lis— Murgida and co-worker a 10 ' 11 presented an 
scheme to reach desired excited states in two-electron 
DQDs by "navigating" in the energy spectrum with elec- 
tric fields - an approach widely used for two- and three- 
level systems Here we focus on a similar quasi-one- 
dimensional setup and show that optical fields can lead to 
almost 100% fidelities of targeted entangled two-electron 
states in tens of picoseconds, thus improving both the 
operation time and the fidelity from earlier works by 1-2 
orders of magnitude. In this context, we extend the previ- 
ous single-particl o 14 ' 15 and two-particlel 6 .^ optimization barrier can be adjusted by a scanning gate microscope 




FIG. 1: (color online) (a) Schematic picture of the pro- 
posed setup. Two quantum dots in a quasi-one-dimensional 
nanowire are separated by tunnel barriers, and the central one 
can be adjusted by a conducting tip. (b) One-dimensional 
model potential for a double quantum dot used in the calcu- 
lations. 



studies into an all-round scheme to coherently control an 
arbitrary transition in a two-electron DQD. Moreover, we 
show that the control procedure validates for the max- 
imization of the entanglement itself, which might have 
further implications for quantum information in solid- 
state devices. 

Our model system is illustrated in Fig. [IJa) . We con- 
sider a quasi-one-dimensional nanowire fabricated, e.g., 
by epitaxial methods* 1 ^ Two QDs are separated by tun- 
nel barriers, and the strength and position of the central 



tip 
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Here we assume the tip to be located in the mid- 



dle, so that the one-dimensional model potential has the 
shape shown in Fig. [IJb). The edges of the barriers are 
softened by Gaussian functions. 

For general applicability of our results, and to explic- 
itly compare with Refs.[l(jand[Il], we consider GaAs ma- 
terial parameters with the effective mass m* = 0.067 m e 
and dielectric constant e = 12.4 eq. We use effective 
atomic units (a.u.) throughout the paper apart from 
the physical model in Fig. [T^b). The energies, lengths, 
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and times scale as E£ = (m* /m )/(e/eo) 2 Eh ~ llmeV, 
a o = i £ / £ o)/( m * / m o) a o ~ 10 nm, andig = fi/E^ sw 55 fs, 
respectively. 

The two-electron problem in one dimension is solved 
numerically exactly by considering electron coordinates 
xi and x 2 and a soft-Coulomb interaction of the form 
Vi n t(xi,x 2 ) = [(xi — x 2 ) 2 + a] -1 / 2 with a softening pa- 
rameter a = 0.25. The Hamiltonian reads H = K + 
V ext (xi : x 2 ) + V lnt (xi, x 2 ) - fie(t), where K is the kinetic- 
energy operator and V ex t(xi, x 2 ) = maji[V(xi),V(x 2 )) 
with V(xi) is the external confinement potential for each 
individual electron shown in Fig. Hlb). It represents a 
double quantum well with a depth 220 meV, width 35 
nm, and interwell barrier 5 nm. In the Hamiltonian 
fx = — {x\,x 2 ) is the dipole operator and e(t) is the op- 
tical field. We point out that the spatial wave functions 
across the diagonal X\ = x 2 are symmetric, since we con- 
sider only antisymmetric two-electron singlet states (so 
that the total wave function is antisymmetric) . We could 
focus on triplets in a similar manner by considering only 
spatially antisymmetric wave functions. A small distor- 
tion potential Vd{x\,x 2 ) = 0.001(xi +x 2 ), corresponding 
to a small static field, is added to the Hamiltonian to dis- 
tinguish degenerate eigenstates as in a real setup, where 
impurities are always present. 

Our task is to optimize e(t) such that the two-electron 
wave function \^(xi,x 2 )} reaches the desired target state 
|$f) at a fixed time T by maximizing the overlap 

|(\I/(xi,X2,T)|<I>f)| 2 - To this end, we apply quantum 
optimal control theory^ (OCT) as explained in Ref. UH. 
This gives the optimized field e(t) through an iterative 
procedure of the time-dependent Schrodinger equation. 
As physical constraints, we set the maximum allowed 
frequency in the field to w max and the fluence F, i.e., 
the integrated intensity (values given below) . We use the 
OCT scheme of Werschnik and Gross^ 2 - and apply the 
octopus code 2 ^ in real space and real time. 

First we consider an excitation where the electrons stay 
in separate dots (we call this a delocalized configuration) 
but a node is created in both dots as shown in the in- 
set of Fig. [21 This corresponds to an excitation from 
the ground state |1) to the fourth state |4) (|2) and |3) 
are localized configurations with both electrons being in 
the left and right well, respectively) The obtained yield 
strongly depends on w max allowed in the optimization, 
and with a pulse length T = 300 (17 ps) we find a strong 
increase slightly above w max = 1.4. This is simply due to 
the resonance at E4 — E\ = 1.3957 (dashed line), which 
is a crucial frequency to be included in the optimization 
to reach extremely high yields. A shorter (T = 100) but 
more intensive pulse lacks this behavior due to the tun- 
neling character of the excitation; the Keldysh parameter 
- used in atoms to assess whether the process is of tun- 
neling or (multi)photon character - is inversely propor- 
tional to the peak intensity.™ However, the shorter pulse 
leads to higher yields than the longer (and less intensive) 
one with filter frequencies smaller than the resonance fre- 
quency. To underline the effectiveness of OCT in Fig. [5] 
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FIG. 2: (color online) Yield as a function of the maximum 
frequency in the optimized pulse to excite the system from 
the ground state to the fourth state. The electrons remain 
delocalized but a node is created in the wells (see the inset). 
Two different pulse lengths and fluences are considered (see 
text). 



it is noteworthy that a non- optimized pulse with the res- 
onant frequency and T = 100 and T = 300, leads, at 
most, to yields ~ 70% and ~ 95%, respectively, whereas 
the optimized pulses, with u; max about 10% above the res- 
onance, lead to respective yields of ~ 95% and ~ 99.95%. 

In Fig. [3] we show our results for an opti- 
mized excitation from the ground state |1) to 
a superposition that consists of the six lowest 
states, i.e., the target state can be written as 
|$ P ) - (i/V6)(|l)+|2)+|3)+|4)+|5)+|6)) 
[Fig. [3{a)]. The pulse length and fluence are fixed to 
T = 600 (= 33 ps) and F = 6 (= 71 meV), respectively, 
and the filter frequency is ujp = 3 (= 8.6 THz). As 
shown in Fig. E^b) we find an extremely high total 
yield 99.6%, and the occupations of the individual 
states converge close to 1/6 (see the inset). The final 
density, consisting of both localized and delocalized 
contributions, is shown in the right inset. 

It was shown in Ref. [ll| that properly adjusted electric 
fields can drive a similar system to a three-state super- 
position in 1000 ps with 96.56% final overlap. In this 
respect, our optical scheme appears to be about 30 times 
faster and an order of magnitude more accurate, even 
though the target system is more complex (six states 
instead of three). The price to pay is the rather com- 
plicated control field in Fig. |3Jc) that, however, satisfies 
realistic constraits and, secondly, can be thoroughly un- 
derstood with the following analysis. It is important to 
note that the optimized field as the one in Fig. |3Jc) cor- 
responds to a local maximum in a search space, i.e., there 
might be other solutions yielding similar, or even better 
outcomes. 

Figure U shows a time- frequency map of the optimized 
pulse in Fig. |3tc) . The amplitudes of the appearing fre- 
quencies (y axis) during the pulse (x axis) are shown 
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FIG. 4: (color online) Time- frequency map of the optimized 
field shown in Fig. [3jc). Amplitudes of the appearing fre- 
quencies as a function of time are shown in colorscale (white: 
small, black: large). Several (de) excitation processes can be 
distinguished (see text). 



FIG. 3: (color online) (a) Excitation process from the ground 
state to an equal superposition of six eigenstates. (b) Occu- 
pations of the states during the optimized process. At the 
end (inset) the states reach ~ 1/6 occupations while a total 
yield of 99.6% is obtained for the superposition (total density 
shown in the right inset), (c) Optimized pulse analyzed in 
Fig- HI 

in colorscale. Time frames of At = T/10 = 60 were 
used in the Fourier transforms. Several distinctive fea- 
tures can be seen in the map, and they can be associated 
with specific (de)excitation processes. First, in the begin- 
ning of the pulse we find pronounced frequencies around 
uj 1.5 and 2.4 corresponding to delocalized excitations, 
i.e., nodes are created in left-right configurations, with 
cji j4 = Ei - E\ = 1.40 (as in Fig. 0), ^4,7 = 1-43, and 
W4.10 = 2.30. Around t = 400 . . . 500 we find similar pro- 
cesses for localized configurations (left-left or right-right) 
with W2,5 = 6 = 1-57 and W5.8 = a_>6.9 = 1.06. In be- 
tween, we find contributions corresponding to transitions 
between delocalized and localized configurations. They 
occur through superpositions of the localized states and 
require frequencies equal to energy gaps between delocal- 
ized and localized states. Finally, close to the end of the 
pulse we find a very large amplitude at uj ~ 0.5 alongside 
several small-scale features. They correspond to several 
deexcitations from high-lying occupations back to the tar- 
get eigenstates as discussed within Fig. [3] Overall, with 
the time-frequency map we are able to understand most 
features of the optimization process despite the compli- 
cated shape of the pulse in the time domain [Fig. EJc)]. 
On the basis of Fig. [2] it can be expected that filtering 
out the dominating frequencies in Fig.Umay significantly 
lower the yield. Quantitative assessment of the sensitiv- 
ity of the yield to different frequencies is the subject of 



future work. 

Finally, we examine how the correlation entropy and 
entanglement are evolved during the controlled excitation 
to a superposition. The correlation entrop y 25 ' 26 is de- 
fined as S — — rifc In nk, where the sum runs over the 
N lowest eigenvalues of the one-body density matrix. 
The maximum value for S is obtained when = 1/N, 
i.e., when all the states are equally occupied. 

The entanglement, on the other hand, is defined 
here in terms of the von Neumann entropy 2 ^ Ej = 
— Tr (pj log 2 pj), where pj = Tij (^1) is the reduced 
density matrix, and the trace is taken over the j:th 
site (1 or 2) of the DQD. The three possible well oc- 
cupations are |q) ■ with q = 0,1,2 being the charge 
of the well. Now the reduced density matrix becomes 
Pj = Po,j (0|j + Pi,j + P2,j \2) j (2\j (occupa- 

tions and \2)j should not be mixed with states |1) 
and 1 2)). Here pij is the probability of the jth site being 
occupied with i electrons. Therefore, the entanglement 
can be written as Ej = ^Zi=oPi,j ^°&2Pi,j- 

In Fig. [5]we plot the correlation entropy (upper panel) 
and the entanglement (lower panel) for two excitations: 
from the ground state to the superposition of the first 
three states, i.e., |$ F ) = (l/%/3) ( |1) + |2) + |3)), and to 
the superposition to the first six states (the same as in 
Figs. El and 0}. The correlation entropy oscillates strongly 
in both cases but eventually reaches the values 1.1112 
and 1.8062, respectively, that are very close to the ideal 
values ln(3) 1.0986 and ln(6) « 1.7918 (dotted lines). 
This means that at t = T the eigenstates in the target 
superpositions are almost equally occupied. The higher 
entropies at t < T are due to the distribution of the 
electron occupations to a large number of states - about 
20 states are occupied by more than 1% in the process. 
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FIG. 5: (color online) Correlation entropy (upper panel) and 
the entanglement (lower panel) in excitation processes from 
the ground state the superposition of the first three states 
(solid lines) and to the superposition of the first six states 
(dashed lines). The latter process is the same as in Figs. [3] 
and |4] The horizontal dotted lines in (a) and (b) mark the 
ideal and maximum values, respectively (see text). 

The entanglement in the lower panel of Fig. [5] ap- 
proaches the theoretical maximum value log 2 (3) ~ 1.5850 
(dotted line). For the three-state and six-state optimiza- 
tion processes we find 1.5764 and 1.5830, respectively. 
This can be associated with the maximization of the en- 
tanglement with 99.5% and 99.9% fidelities, that are com- 
parable to the respective target-state overlaps of 99.5% 
and 99.6%. Thus, through the proposed setup we are 
able to find maximally entangled two-electron states in 
a targeted fashion. We also note that the entanglement 
could be explicitly measured in the DQD by detecting 
the charge^ on one of the wells leaving the other well 
unaffected. This procedure would directly correspond to 
tracing out one site of the DQD and reducing the density 
matrix of the whole system to pj = Trj (\&|). 

We point out three aspects in the practical realization 
of the proposed control scheme. First, the actual con- 
fining potential and thus the Hamiltonian could be de- 
termined through a supporting experiment, e.g., by mea- 
suring the single-particle spectrum and constructing V ex t 
through inversion.— The optimization could then be per- 



formed within this specific V^xt . Secondly, it is important 
to note that the physical constraints of the optical pulses 
required to induce the proposed transitions are within 
reach of the present experimental capabilities. Terahertz 
frequencies can be routinely produced with, e.g., quan- 
tum cascade lasers and they can be shaped even in the 
femtosecond time scale.— For example, in a recent ex- 
periment^ femtosecond pulses were created in a synthe- 
sized fashion by separately controlling the infrared, vis- 
ible, and ultraviolet contributions, including their chirp, 
carrier envelope phase, delay, and energy (beam size). 
Therefore, coherent control with THz optics is likely to 
become a feasible scheme in solid-state quantum infor- 
mation. On the other hand, voltage control with gates is 
so far strictly limited to GHz time scales. Finally, as the 
third critical aspect we mention coherence that highly 
depends on the external conditions in a particular de- 
vice. We believe that the time scales presented here are 
coherently reachable, although we cannot rule out the 
possiblity of decoherence due to hyperfine effects or/and 
interactions with optical and acoustic phonons. 

To summarize, we have shown through numerically 
exact calculations that coupling a two-electron double 
quantum dot to an optimized optical source enables co- 
herent excitations to entangled states with an extremely 
high fidelity In the proposed setup the pulses produced 
by quantum optimal control theory can be exposed to 
realistic constraints without losing the accuracy. Time- 
frequency analysis gives invaluable insight on the physical 
processes in the evolution towards the desired superpo- 
sition. We show that entanglement of delocalized and 
localized electron configurations can be reached up to 
99.9% of the theoretical maximum. We hope that the 
proposed scheme motivates further experiments in the 
initialization, control, and readout of electric charge in 
solid-state devices. 
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